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Abstract 

We study the signal of DM and photon associated production induced by the dimension six vector 
and axial-vector operators at the LHC, including the QCD next-to-leading-order (NLO) corrections. 
We find that the QCD NLO corrections reduce the dependence of the total cross sections on the 



factorization and renormalization scales, and the K-factors increase with the increasing of the DM 
mass, which can be as large as about 1.3 for both the vector and axial-vector operators. Using our 
QCD NLO results, we improve the constraints on the new physics scale from the results of recent 
CMS experiment. Moreover, we show the Monte Carlo simulation results for detecting the 7 + Ifj, 
signal at the QCD NLO level, and present the integrated luminosity needed for a So" discovery at 
the LHC with ^/S = 14 TeV. If the signal is not observed, the lower limit on the new physics scale 
can be set. 

PACS numbers: 12.38.Bx, 14.65. Jk, 14.70.Bh, 95.35. +d 



* Electronic address: csli@pku.edu.cn 



1 



I. INTRODUCTION 



The dark matter (DM) attracts a lot of attention in the fields of both cosmology and 

n n 

particle physics [1|, |2| . The astrophysical observations have provided strong evidence for 
the existence of DM [3| . The relic abundance of DM favors the weakly interacting massive 
particles and this kind of DM has been extensively studied. Since there are no DM candidates 
in the Standard Model (SM), any discovery of the DM signal implies new physics (NP) 
beyond the SM. 

There are mainly three kinds of methods to detect the DM. First, we can detect DM 
annihilation signals from galaxy clusters, such as PAMELA Q, ATIC y, HESS y and 



Fermi-LAT 0, 18| . Second, we can search for the signal of DM underground when the DM 



^articles interact with nuclei, such as the experiments of XENON jo], CDMS 



10|, DAMA 



ll| and CoGeNT 12|. Third, we can study the properties of DM at hadron colliders if 



DM exists and has interactions with the SM particles [13|, [l^ . Compared to the direct and 
indirect experiments, the hadron colliders have impressive advantages that the measurements 
are not sensitive to the uncertainties related to the galactic distributions, DM velocities, etc. 
There have been a lot of studies to search for DM at the LHC in a series of DM models 



15 



32|. 



At the LHC, the DM manifests itself as missing energy when produced. However, we 
can probe the DM through the visible particles, which are associated produced, such as a 
photon or a jet 33|, |34]. In this work, we only consider the DM and photon associated 
production at the LHC, since this signal is clear and suffers from less backgrounds from the 
SM. Since the fundamental theory of the DM is unknown, it is useful to investigate the DM 
in a model-independent way by using effective field theory. 

Recently, the CMS collaboration has searched for NP in the 7 + Ifj. final state, and 
set the 90% confidence level (CL) lower limits on the NP scale for vector and axial-vector 
operators 13|]. However, the analysis there is based on the leading order (LO) results, which 
suffer from large uncertainties due to the choice of renormalization and factorization scales. 



Following our previous works 35|, |36(|, in this paper, we study the signal of DM and photon 



associated production induced by the dimension six vector and axial-vector operators at the 
LHC, including QCD next-to-leading-order (NLO) corrections. Using our NLO results, we 
improve the constraints on the NP scale from the results of recent CMS experiment. 
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The paper is organized as follows. In Sec. HIl we describe the dimension six vector and 
axial-vector operators describing the interactions between DM and the SM particles. In 
Sec. Illlt we calculate the relic abundance induced by these effective operators and find the 
allowed regions for the mass and coupling of the DM. In Sec. IIV[ we present the details 
of the QCD NLO corrections to the associated production of the DM and photon at the 
LHC In Sec. |Vl we discuss the scale dependence and show the QCD NLO cross sections for 
different DM masses numerically. After that we improve the experimental constraints given 
by CMS collaboration. In Sec. IVIt we calculate the backgrounds in the SM and analyze the 
discovery potential at the LHC with = 14 TeV. Finally, we conclude in Sec. IVIII 



II. EFFECTIVE OPERATOR 



We consider the dimension six vector and axial-vector operators 



Ov 

Oa 



K 

A2 

K 

A2 



(1) 



which are also studied in Refs. 



25 



33 



37 



38|. The NP scale A can be regarded as 



the remnant of integrating the massive propagator between the DM and SM particles. We 
assume that the Dirac fermion x is a DM candidate, and a singlet under the SM gauge group 
SU{3)c X SU{2)i X U{1)y- The DM x can only interact with the quarks by these operators. 



In Ref. 



13| . the constraints on the NP scale A are given by the CMS collaboration through 



the process of DM and photon associated production at LO. In this paper, we will perform 
the QCD NLO corrections to these processes, whose effects are important for research at the 
LHC, and improve the limits on the NP scale. For simplicity, we only present the details of 
the calculations for the vector operator. For the axial operator, we only give some necessary 
analytical results and the numerical results. The Feynman diagrams are the same for both 
vector operators and axial-vector operators except for the vertices. 



III. RELIC ABUNDANCE 



Before discussing the signal of the DM at the LHC, we first consider the relic abundance 
which is a precise observable in cosmology. The DM we choose to study can contribute to 
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(c) 

FIG. 1: Feynman diagrams for the DM annihilation process. Labels (a), (b) and (c) denote the 
LO, virtual correction and real corrections, respectively. 

the relic abundance of cold dark matter (CDM). Therefore, the relic abundance imposes 
constraints on the DM. Here we use the experimental results for the relic abundance in Ref. 



39| 



VLcDMh^ = 0.1123 ± 0.0035, 



(2) 



where ^Icdm is the CDM energy density of the universe normalized by the critical density 
and h = 0.710 ± 0.025 is the scaled Hubble parameter. 

The relic abundance can be calculated from the total annihilation cross section of DM, for 
which the Feynman diagrams are shown in Fig. [H First, we give the LO total annihilation 
cross section for the vector operator 

' s s + 2m^ 



'B,V 



f- 



(3) 



A4 V s - 4m2 127r ' 

where s is the usual Mandelstam variable. A^^^ and Nj are the numbers of color and flavor 
of quarks, respectively, m is the mass of the DM. This LO cross section is consistent with 
the result in Ref. js?]. 

The QCD NLO corrections to the total annihilation cross section include two parts, i.e., 
the one-loop virtual corrections and real gluon emission corrections. We use n = 4 — 2e 
dimensional regularization to regulate the divergences in the calculation. The results of 
virtual corrections are 



a. 



vy 



a 



By 



OsCf 
271 



(4) 
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where = (47r/i^/s)^pp3|^ and Cp = 4/3. The resuhs of real corrections are 

'asCF\ „ /2 3 19^ 



+ - + ^ . (5) 



'ry - -By \ 2n J^'\e^ ' e ' 2 
Combining the two parts, we get the QCD NLO total annihilation cross section 

(^NLoy = ^ (^Byj W 

in which K"'"' is the fC-factor, generally defined as ctnlo/clo; of the DM annihilation cross 
section 

i^«- = l + lfi. (7) 
vr 

The DM is moving at nonrelativistic velocities {v ^ 1) when freezing out. We define v 
as the relative velocity between the DM. Thus we can expand 

(^NLoyV = a + bv'', (8) 

where 

b = K'^^N.Nf'^^-^. (9) 

^ A4 247r ^ ' 



We perform the calculation of relic abundance by using the method in Ref. 40|. The 
freeze-out epoch xj = m/Tf is given by 



X/ = In [0.038ampim {9/9^^)] - ^ In {in [omSampim ((7/(7^)] } 
+ In < 1 + - 



a In 



0.038ampim [9/9 



1/2 



(10) 



where mpi = 1.22 x 10^^ GeV is the Planck mass. 9 and is the number of relativistic 
degrees of freedom and effective number of relativistic degrees of freedom at the freeze-out 
temperature Tf, respectively. The result for the relic abundance is 

1.07 X 10^ GeY-'x, 
{9*3/9* )fnpia 

For most of the history of the universe all particle species had a common temperature, and 
9^:S can be replaced by 9^:. Requiring the DM relic abundance is in the 2a region around the 
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FIG. 2: Relation between the mass of DM and the NP scale for the vector operator. The 
relic abundance is required be to in the 2a region around the observed central value. The lower 
green band is the LO result. The upper red band is the NLO result. In this figure, we choose 
0.118 and Nf = 5. 



K = l,as 



observed central value, the NP scale of the effective operator is determined by the mass of 
the DM. From Fig. [21 we see that the NLO corrections increase the lower limits on the NP 
scale slightly. Since we do not hold the point of view that the relic abundance is determined 
only by this one kind of DM, the regions below the red band are all allowed. 

For the case of axial-vector operator, we follow the same process and give only the main 
results. The born cross section is 

s s — 4m^ 



N.N, 



K 



(12) 



After taking the QCD NLO corrections to the process induced by the axial-vector operator, 
we get 



(13) 



in which 



a 
b 



(14) 
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FIG. 3: Relation between the mass of DM and the NP scale for the axial-vector operator. The 
relic abundance is required be to in the 2a region around the observed central value. The lower 
green band is the LO result. The upper red band is the NLO result. In this figure, we choose 
K = = 0.118 and Nj = 5. 

Then we get the constraints from the relic abundance as shown in Fig. [3] by using the 
following formulae 

X/ = In [Q.mhmpim [g/g]'^)] - ^ In {in [0.4566mp;m [g/g]''')] ] , (15) 



1.07 X lO'^ GeV^'x 



-1^2 
/ 



1/2 • (16) 

^{g*s/g* )mpib 

From Fig. |3l we get nearly the same results that the NLO corrections increase the lower 
limits on the NP scale slightly, and the regions below the red band are all allowed. 



IV. QCD NLO CORRECTIONS TO DM AND PHOTON ASSOCIATED PRO- 
DUCTION 

From the effective operators in Eq. ([T]), the DM can be pair produced via annihilation of 
quark pairs. However, because DM can not decay into SM particles, such processes give just 
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FIG. 4: LO Feynman diagrams of the DM and photon associated production. 



missing energy and no observable signals. As a result, we have to consider the process of DM 
associated production with a photon or a jet 



25 



33 



and use the initial radiation to tag 



the DM. Though the cross section of associated production with a photon is less than that 
of associated production with a jet, the signal of this process is characteristic and involves 
less backgrounds. 

Different from the scalar operator in our previous work [ssl, the quark sector and the DM 
sector can not factorize for the vector and axial-vector operators that we consider in this 
paper. This leads to more complicated analytical expressions in our calculation. We only 
show the analytical results for the vector operator in detail and give the numerical results 
for both of the vector and axial-vector operators. 



A. LO calculation 



First of all, we calculate the LO cross section of the following process 

q{pi) + q{p2) x{P3) + X{P4) + l{P5)- (17) 

The LO Feynman diagrams are shown in Fig. |H Then the LO partonic cross section is 

1 



2s 



12 



dTs\Mi 



(18) 



in which Fs is the three-particle final states phase space. We define Sij = {pi tij = 

{Pi ~ PjY ^'^d ^ — e^/47r. The spin and color summed and averaged Born matrix element 
squared is 

.2 



\Mb? 



rran 



2| \^V{A),2 
I ' 



(19) 



in which |A^^P is expressed as 

1 



\M^o\ = -— 16(2m^(4si2 + tl5 + t25)+m"(-Sl2(2S35 
^15^25 

+2S45 + 3tl3 + 3tl4 - 4ti5 + 3t23 + 3t24 " 4t25) 
"■S 45^15 — S45i(:25 ^ S35(tl5 + ^25) + 

4(^12)^ + 2(ti5)^ + 2(t25)^ ~ 2ti3ti5 — 2ti4ti5 + ti5t23 + ^15^24 + 
^13*^25 + ^14^25 — 2^23^25 ~ 2^24*^25) + ■545ti3ti5 + 
535^14^15 + 512(545(^13 + ^23) + 

535(^14 + ^24) + 2(ti4t23 + ^13^24)) + ■545t23^25 + 535^24^25 + 
^141^15^23 + i^l3^15'^24 ~ '^'tlb't 23^24, ~ '2'tl3'tu't25 + 

^141^23^25 + 1^13^241^25)) (20) 

and Qi{i = 1, 5) are the electric charge of the quarks. For the axial-vector operator, 

|_y^A|2 ^ |_^y|2 _ 64m^ (2^12 (tis + %) + 2 (^12) ^ + (^15) ^ + (%) ^) ^^i) 

^15^25 

After convoluting with the parton distribution functions (PDFs) Gg(q)(x), we obtain the 
LO cross section 

(Tb = y dxidx2[Gg/p{xi)Gq/p{x2) + (xi ^ X2)]aB- (22) 
B. QCD NLO corrections 

The QCD NLO corrections consist of real gluon radiation, quark or antiquark emission 
and one-loop virtual gluon effects. We use n = 4 — 2e dimensional regularization to regulate 
both the ultraviolet (UV) and the infrared (IR) divergences in our calculations. 

The virtual gluon corrections involve the self-energy, vertex and box diagrams, as shown 
in Fig. [5l We use on-shell renormalization scheme for external fields. The UV divergences 
are exactly canceled between the loop diagrams and the counterterms of the external fields. 
For the vector operator in our case, there is no renormalization for the new coupling in this 
renormalization scheme, same as the Drell-Yan process. We denote the bare and renormal- 
ized quark wave functions by ipgo and ipq, respectively. The renormalization constant 6Zg is 
then defined by 

lljgO = {l+SZgy/'^g. (23) 




FIG. 5: Feynman diagrams for one-loop virtual corrections. 



Calculating the quark self-energy diagram, we obtain the explicit expression for 6Zq: 



(24) 



Here 1/e and l/ej/y represent IR and UV divergences, respectively. 

After renormalization, the UV divergences in the virtual corrections are removed, leaving 
the IR divergences and the finite terms. The final virtual gluon corrections to the partonic 
cross section are 

(25) 



dT32Re{M*BM^). 



The IR divergent part of A^^ is given by 



B, 



(26) 



where = r(l + e)[(47r//|)/si2]' and 



Al 



-3Cf. 



(27) 



Another possible contribution at hadron collider arises from the process 



9 + 9 ^ X + X + l, 



(28) 



which is illustrated in Fig. |6l The effects of these loop-induced gg diagrams are one order 
of as higher than that of qq diagrams. The triangle diagrams in Figj6] vanish because of the 
color structure. The box diagrams do not vanish, but the gluon PDF is strongly suppressed 
if we require the final states have large transverse momentum. Thus, we assume their 
contribution is small and do not discuss this process in the rest of this paper. 
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FIG. 6: Sample Feynman diagrams for gluon-gluon initial states contributions. 




fWW 



^AAA/ 



fVWV 



AAAA/ 




FIG. 7: Feynman diagrams for a real gluon emission. 
The Feynman diagrams for the real gluon radiation process 

q{pi) + q{p2) xiPs) + X{P4) + 7(^5) + 9{Pg) 



(29) 



are shown in Fig. [71 Soft and colhnear divergences appear when we perform the final state 
phase integrations. To cancel the IR singularities, we iise the two cutoff phase space slicing 
method to integrate the singular regions analytically [41]. Explicitly, we use the soft cutoff 
parameter Sg to define the soft regions and the collinear cutoff parameter 6c to define the hard 
collinear regions. The soft regions are just the phase space where the real radiated gluon's 
energy Eq < ds^/syi/'^- The collinear regions are defined by Itjgl < 5cSi2 with i = 1,2. 
Therefore, the real corrections are divided into three parts, i.e.. 



ddr = da^ + da^^ + da. 



■HC 



(30) 



Here, and a^'~^ represent the partonic cross section for the soft regions and hard collinear 
regions, respectively. The hard non-coUinear part a^*" is finite and can be computed nu- 
merically. 

In the soft regions, we use the Eikonal approximation, so that the cross section can be 
factorized as 

daf. = {A7rasfJ,R)daB J dS^eik, (31) 



11 



where the phase space integration factor dS over the soft region is expressed as 

'•Ss^/Sl2/2 



dS 



And the Eikonal factor ^eik is the amphtude squared in the soft hmit apart from |A^sP, 
expressed as 

^eik^CF-^. (33) 

After integration over the soft gluon phase space, we get 

da^ = das^C, + T + ' ^^^^ 

where 

AI^2Cf, Af^-ACplnS,, ^ (^A\n' 5^ - (35) 

In hard colhncar regions, the momentum of the gluon emitted from initial partons become 
coUinear to the beam line. In this limit, the four-body matrix elements are approximated 
as follows: 



|A4^|2 {A7rasiiR)\MB\ 



(36) 



ZtiQ Zt26 

in which z represents the fraction of initial partons' momentum carried by q{q). Pij{z,e) 
arc the unregulated splitting functions in n-dimensions which can be related to the usual 
Altarelli-Parisi splitting kernels as Pij{z, e) = Pij{z) + eP-j{z). In our case, 

PU^) = Cf\^, P',,{z) = -Cf{1-z). (37) 
Because of charge conjugation invariance we have 

Pqq = Pqq, Pqq = Pqq- (38) 

At the same limit, the four-body phase space can be written as 

dTi\coii = rfr3(s^2 = ^^12) 2g^2r("j _ ^^ dzdtiQ[-{l - z)tiQ]''. (39) 
Therefore, we obtain 

dcrf^ = d&B^C, (^-^^ S~'[Pqq{z,e)Gq/p{Xi/z)Gq/p{x2) 

dz fl - z\~^ 

+ Pqq{z, e)Gg/p{xi)Gq/p{x2/z) + {xi -H- X2)]— ( j dxidx2. (40) 

12 



To factorize the collinear singularity into the PDFs, we use scale dependent PDFs in the 
MS convention: 



Gb/p{x,HF) = Gh/p{x) + 



a. r(i 



/ 47r/i' 



^ dz 

-PU^)GaUx/z). (41) 



_27rr(l -2e) V 1^1 

Now, we replace Gq(^q)/p in the LO hadronic cross section (l22l) and combine the result with 
the hard collinear contribution in Eq. ( l40l) . The resulting 0{as) expression for the initial 
state collinear contribution is 



rfcr™" = daB^G^\Gq/p{Xi,^lF)Gq/p{x2,llF) + Gq/p{Xl,^F)Gq/p{x2,^lF) 



a. 



Af{a — )■ ag) 



+ Al''{a ag) Gg/p{xi, fiF)Gg/p{x2, ^if) 



a=q,q 



+ (Xi X2)^dXidX2- 



(42) 



with 



AT{q^qg) = Cf(21n5, + 3/2), 
A'o\q-^qg) = ATiq ^ qg) In T-^). 

The G functions are given by 



Gb/pix,lJ.F) = 



l-5sSab 



dy ~ 

— Ga/p{x/y,HF)Pba{y) 

y 



(43) 



(44) 



with 



(45) 

A complete real correction includes also the (anti) quark emitted processes, as shown in 
Fig. m such as 

9{V\) + qlq{V2) xiPs) + X{P4) + 7(^5) + q/qiPe)- (46) 

Their contributions can be obtained from the results of processes in fl29p by crossing sym- 
metry. Note that, when separating the singular phase space regions for these processes, one 
only needs to deal with the collinear divergences which can be totally absorbed into the 
redefinition of the PDFs in Eq. (14T|) . 

Finally, the NLO cross section for the process pp — ?■ XXI is 

- J dxidx2{ [G,/p{xi,fiF)G^/p{x2, fip) + {xi ^ X2)] {as + + + af^)} + (t^°" 
dxidx2[Gg/p{xi, i^F)Ga/p{x2, fip) + {xi <r^ X2)](Jr'^' {ga XXia), (47) 

a=q,q 



a 
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FIG. 8: Feynman diagrams for a quark emission. The Feynman diagrams for antiquark emission 
can be obtained by charge conjugation. 

where C in dr'"{ga — )■ XXI'^) suggests that the phase space integration is performed in the 
noncoUinear regions. We have checked that 

+ = 0, A\ + Al + 2Af{q -> qg) = 0. (48) 

Therefore there are no singularities left any more and we can perform numerical integration 
now. 



V. NUMERICAL RESULTS 

In this section, we give the numerical results for the cross sections of DM and photon 
associated production at the LHC In numerical calculations, we choose the CTEQ6L1 
(CTEQ6M) PDF sets [42] and the corresponding strong coupling as for the LO (NLO) 
calculations. The default factorization and renormalization scales, fip and fiR, are set as 
2m. Recently, the observations of the gamma-ray in Fermi-LAT gives the hints of 130 
GeV DM Thus, we choose the input parameters (m,A) = (130 GeV, 500 GeV) and 
K = 1 unless otherwise specified, which are allowed by the relic abundance constraint. The 
kinematic cuts 

pI > 100 GeV, 
\V^\ < 2.4, 

p^''' > 100 GeV, (49) 
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FIG. 9: Dependence of the NLO cross sections for the DM and photon associated production at 
the LHC on the soft cutoff 6s with 6c = 6s/ 50 for the axial- vector operator. 

are applied in our numerical calculations. Here p^*^* is the missing transverse momentum, 
defined as 

{nZ, no iets in the final states, 
, with jets in the final states, 

where p^^ is the transverse momentum of the system of the DM. Jets are defined by the 
requirements pi^^ > 20 GeV and \ri^^^\ < 2.5. In order to avoid QED collinear divergences. 



we also require the photon to be isolated by the prescription 



43| 



where R = ^/Acp^ + Ar]'^ and Rq = 0.4. 

Since we use the two cutoff phase space slicing method in this calculation, it is necessary 
to show that the total cross section is independent of the cutoff parameter. In Fig. [9l the 
contribution of three-body final states includes the Born cross section and one-loop virtual 
corrections, soft and collinear limits of the cross section of four-body final states. The 
contribution of four-body final states consists in the cross section of four-body final states 
with the singular regions of the phase space sliced. From Fig. [9l we see that our final result 
is nearly independent of the cutoff parameter, since the total cross section varies less than 
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FIG. 10: Dependence of the LO (NLO) cross sections for the DM and photon associated production 
at the LHC on the factorization scale hf and renormalization scale nn for the vector operator. 




1.2 1.4 1.6 1.8 2 0.6 0.8 1 1.2 

|i/2m [i/2m 



FIG. 11: Dependence of the LO (NLO) cross sections for the DM and photon associated production 
at the LHC on the factorization scale /ip and renormalization scale /ir for the axial- vector operator. 

5% when the parameter 6s changes from 10^^ to 0.005. 

In Fig. [To] and Fig. [TTl we show the dependence of the LO (NLO) cross sections for 
the DM and photon associated production at the LHC on the factorization scale fip and 
renormalization scale /i/j. It can be seen that the dependence of the NLO cross section on 
the factorization scale fip and renormalization scale is significantly reduced, compared 
to the LO cross section. This makes the theoretical prediction much more reliable. 

In Fig. [121 we show the DM mass dependence of the LO and NLO cross sections for the 
vector operator with \/S = 14 TeV. When the DM mass varies from 130 GeV to 200 GeV, 
the QCD NLO corrections are modest. For the DM mass in the range from 300 GeV to 1000 
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GeV, the QCD NLO corrections generally improve the cross section and are more significant 
for larger DM mass. For example, the QCD NLO corrections increase the cross sections by 
about 19% for m = 1000 GeV. Thus, it is necessary to consider the NLO corrections to the 
process of DM production at hadron colliders. We also show the mass dependence of the 
K-factors for the axial-vector operator in Fig. [T31 which are similar to the case of vector 
operator. In order to compare with the experimental results of CMS, we use the same cuts 
and center-of-mass energy as in |l^, and can improve the lower limits on the NP scale in 
the results of CMS [13], using our K-factors at the LHC with = 7 TeV. Here, we show 
the improved limits on the NP scale A in Table Hand Table HTl for the vector and axial- vector 
operators, respectively. 

TABLE I: Sample results of the 90% CL lower limits on the NP scale A for the vector operator. 
The LO results are given in the CMS analysis Q. The K-factors at the LHC with VS = 7 TeV 
for different DM masses are also shown. 

m [GeV] A [GeV](LO) [13] A [GeV] (NLO) K-factor@7 TeV 

200 549 564 1.11 

500 442 463 1.20 

1000 246 263 1.31 



TABLE 11: Sample results of the 90% CL lower limits on the NP scale A for the axial-vector 
operator. The LO results are given in the CMS analysis [l3]. The K-factors at the LHC with 
^/S = 7 TeV for different DM masses are also shown. 



m [GeV] A [GeV](LO) [13] A [GeV] (NLO) K-factor@7 TeV 
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VI. BACKGROUNDS AND DISCOVERY POTENTIAL 



In this section, we present the Monte Carlo simulation results for detecting the 7 -|- I^j. 
signal at the QCD NLO level. The main irreducible SM backgrounds are the pp — )• 
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FIG. 12: Dependence of the LO and NLO cross sections on the DM mass for the vector operator 
with = 14 TeV. The K-factors are also shown. 

Y2) + 7 and pp — ?► vv) + j when the jet is misidentified as a photon. We use MCFM 
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47| to calculate the backgrounds at the NLO level. Since we focus on the experiments 



at the LHC, we set the missidentified probability P^jj = 10~^, as pointed in the Ref. 48 1. 

Figure [HI and Figure fTSl show the differential cross sections as functions of p^ and p™^*, 
respectively, for the signal and backgrounds at the NLO level. We can see that the distribu- 
tions of the backgrounds decrease faster than that of the signal as the transverse momentum 
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FIG. 13: Dependence of the LO and NLO cross sections on the DM mass for the axial- vector 
operator with ^/S = 14 TeV. The K-factors are also shown. 

of the photon increases. Thus, the ratio of signal and background will increase if we set a 
larger px cut. 

Figure [16] shows the differential cross sections as a function of r]'' for the signal and the 
backgrounds at the NLO level. We see that the distribution of the signal is more concentrated 
in the central region than the backgrounds. These distributions give some clues to suppress 
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FIG. 14: Dependence of the differential cross section on pj.. 




FIG. 15: Dependence of the differential cross section on p™***. 

the backgrounds more efficiently at the LHC 

Figure [17] presents the integrated luminosity needed to discover the signal at a 5a 
{S/y/S + B = 5) level at the LHC with = 14 TeV. We find that the needed inte- 
grated luminosity grows with the increasing of the NP scale, and depends more strongly on 
the DM mass for larger NP scale. In particular, for A = 1000 GeV and m = 200 GeV, the 
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FIG. 16: Dependence of the differential cross section on rf . 

needed integrated luminosity is 12 fb^^ at the 14 TeV LHC. Fig. [18] shows the results for 
the axial-vector operator. 

In Fig. [ini we present the limits of the NP scale for 3cr {S/\/B = 3) and Sa {S / \/B = 5) 
exclusions at the LHC with = 14 TeV, assuming m = 130 GeV. We see that the NP 
scale is constrained to be larger than 1200 GeV if the LHC with \^ = 14 TeV does not 
detect this signal after collecting an integrated luminosity of 10 fb~^. Figure |20] gives the 
results for the axial-vector operator. 



VII. CONCLUSION 



We have investigated the signal of DM and photon associated production induced by 
the dimension six vector and axial-vector operators at the LHC, including the QCD NLO 
corrections. We find that the QCD NLO corrections significantly reduce the dependence of 
the total cross sections on the factorization and renormalization scales, and the QCD NLO 
corrections are more significant for larger DM mass for both the vector and axial-vector 
operators. Using our NLO results, we improve the constraints on the NP scale from the 
results of recent CMS experiment. Moreover, we calculate the dominant SM backgrounds at 
the NLO level, and show the differential cross sections of both the signal and backgrounds 
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FIG. 17: The integrated luminosity needed for a 5a discovery as a function of the DM mass 
at the LHC with = 14 TeV for vector-operator. We choose the cuts > 300 GeV and 
p-miss ^ QgY above analysis of Fig. [U] and Fig. [15l 
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FIG. 18: The integrated luminosity needed for a 5a discovery as a function of the DM mass at 
the LHC with -v/S* = 14 TeV for axial-vector operator. We choose the cuts pj, > 300 GeV and 
pTniss ^ QgY (jyg ii^Q above analysis of Fig. [U] and Fig. [151 
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FIG. 19: The limits of the NP scale for 3a and 5a exclusions at the LHC with -v/S* = 14 TeV, 
assuming m = 130 GeV for the vector operator. We choose the cuts > 300 GeV and p™*** > 
300 GeV due to the above analysis of Fig. [H] and Fig. [151 

as functions of p^, p™"^*^ and r]'^. The character of these distributions can help to select the 
events in the experiments. Finally, we show the potential to discover the DM at the 14 TeV 
LHC, and provide the exclusion limits on the NP scale if this signal is not observed. 
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FIG. 20: The limits of the NP scale for 3a and 5a exclusions at the LHC with -v/S* = 14 TeV, 
assuming m = 130 GeV for the axial-vector operator. We choose the cuts > 300 GeV and 
pmiss y above analysis of Fig. [H] and Fig. [15l 
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